EXCITATION OF NONLINEAR TWO-DIMENSIONAL WAKE 
WAVES IN RADIALLY-NONUNIFORM PLASMA 

Arsen G. Khachatryan 

Yerevan Physics Institute, Alikhanian Brothers Street 2, 
Yerevan 375036, Armenia 
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The progress in the technology of ultrahigh intensity lasers and high cur- 
rent relativistic charged bunch sources permits the use of laser pulses [1] 
or charged bunches [2] for excitation of strong plasma waves. The excited 
plasma waves can be used both for acceleration of charged particles and 
focusing of bunches to get high luminosity in the linear colliders [2]. At 
present the plasma-based accelerator concepts are actively developed both 
theoretically and experimentally (see the overview in Ref. [3] and references 
therein). 

The amplitude of longitudinal electric field E m3/X in relativistic wake- waves 
excited in cold plasma is limited by the relativistic wave-breaking field [4] 
E rd = [2(7 - l)} 1/2 E WB /(3, where 7 = (1 - /3 2 )~ 1/2 is a relativistic factor, 
P — Vph/c is a dimensionless phase velocity of the wave, E WB = m e u} pe v p h/e 
{E WB [V/cm\ « 0.96n p [cm 3 ]) is the conventional nonrelativistic wave- 
breaking field, ojp e = (47m p e 2 /m e ) 1//2 is the electron plasma frequency, n p is 
the equilibrium density of plasma electrons, m e and e are the mass and abso- 
lute value of the electron charge. The acceleration rate in the wake fields can 
reach tens of GeV/m, that much exceeds the rates reached in conventional 
accelerators. 

The linear wake- field theory is valid when E max -C Ewb- In case of wake 
wave excitation by relativistic charged particle bunch [Plasma Wakefield Ac- 
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celerator (PWFA)] this corresponds to the condition a = \q\ n b /en p <C 1 [2], 
where q is the charge of bunch particles, is their concentration. In the 
scheme of wake wave excitation by a short laser pulse (of the length compa- 
rable with plasma wavelength A p = 27rv p h/ou pe ; Laser Wakefield Accelerator 
(LWFA)), the linear theory is valid when a 2 = e 2 El/{m 2 c 2 u)'l) 1 [5], 
where o;o uj pe and E are the frequency and amplitude of laser radiation 
respectively. The phase velocity of the wake wave is equal to the bunch ve- 
locity Vb in PWFA and to the laser pulse group velocity v g ~ c(l — u 2 e /2ijl) 
( that corresponds to 7 = 7^ ~ cuo/u pe ^> 1) in LWFA. 

The one-dimensional nonlinear wake waves excited by wide drivers (when 
kpTd ^> 1, where k p = u pe /v p h is the wavenumber and r d is the radius of 
the driver) are studied in sufficient detail both for PWFA [6] and for LWFA 
[7]. These studies testify to the feasibility of excitation of strong nonlinear 
plasma waves with the amplitude of up to E re i by bunches with a > 0.5 and 
laser pulses with a 2 > 1. The other important results of the one-dimensional 
nonlinear theory are the steepening of the wake wave and the increase of 
wavelength with amplitude. The wave with the amplitude E max m E re i has 
the wavelength nearly 7 1 / 2 times as large as the linear wavelength \ p . 

In reality, the transverse sizes of the drivers used are ordinarily compa- 
rable or less than their longitudinal size. The allowance for finite transverse 
sizes of the drivers and, accordingly, the transverse motion of plasma elec- 
trons complicate the treatment of the problem in the nonlinear regime. In 
the general case the analytical solution of this regime seems impossible and 
here the use of numerical methods are usually required. The numerical inves- 
tigation of nonlinear effects in two-dimensional (axially-symmetrical) wake 
waves exited in uniform plasma has shown that in the nonlinear regime the 
wavelength becomes dependent on the transverse (relative to the driver prop- 
agation direction) coordinate [3,8-10]. The change of a nonlinear wavelength 
in the transverse direction is due to the dependence of wavelength on the 
amplitude, that, in its turn, is varied in the radial direction owing to finite 
cross section of the driver. This leads to a curvature of the phase front of 
the nonlinear wave [8-10], to steepening and "oscillations" of the field in the 
transverse direction [10,11] and eventually to the development of turbulence. 
From the viewpoint of acceleration and focusing of charged bunches in the 
wave, the curvature of the nonlinear wave front is undesired as the quality 
(emittance, monochromaticity) of the driven bunch worsens. In the present 
work we show that by means of wake wave excitation in plasma, the den- 
sity of which is properly varied in the transverse direction, one can eliminate 
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the nonlinear change of wavelength in the transverse direction and the re- 
lated curvature of phase front. The plasma of this kind may be produced 
by charged beams [12] or laser pulses [13] passing through a neutral gas or a 
partially ionized uniform plasma due to an additional ionization. 

Equations for non-zero components of plasma electrons momentum and 
electromagnetic field describing the steady nonlinear wake-fields in radially- 
nonuniform plasma can be obtained by simple generalization of equations for 
uniform plasma [10,11,14,15]: 

V ± H e + /3— ^ + f3 z N e + (3a = 0, (4) 
dz dz dr 

dE 

N e = N p (r) -a- V ± E r - (6) 

As usual, Eqs. (1) and (2) were derived taking into account the conservation 
of generalized momentum /5 2 H — rotP = 0, or in our case 

In Eqs. (1) - (6) j, = (1 + P* + Pi + a 2 /2)V2 = P z r / 7 , and N e = 

n e /n p (r = 0) are respectively a relativistic factor, dimensionless components 
of velocity and dimensionless density of plasma electrons, N p = n p (r)/n p (0), 
P = v p h/c, z = k p (r = 0)(Z — v p ht) t V ' ±_ = d/dr + 1/r. Also the following 
dimensionless variables have been used: the space variables are normalized 
on X p (r = 0)/2n = l/k p (r = 0), the momenta and velocities - respectively 
on m e c and the velocity of light and the strengths of electric and magnetic 
fields - on the nonrelativistic wave-breaking field at the axis Eywir = 0) = 
m e uj pe (r = 0)v p h/e. The field of forces acting on relativistic electrons in the 
excited field is F(-eE z , -e(E r - /3H e ),0). In PWFA a ^ 0, a 2 = and in 
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LWFA a = 0, a 2 ^ (in Eqs. (l)-(6) the linear polarization of the laser 
pulse field is assumed; for the circular polarization the value of a should be 
multiplied by the factor 2 1 / 2 ). 

We have solved Eqs. (l)-(6) numerically choosing the Gaussian profile of 
the driver both in longitudinal and transverse directions: 

A(z,r) = A exp[-(z - z ) 2 /a 2 z ]exp(-r 2 /a 2 r ), (8) 

where A(z, r) stands for a = nb(z, r)/n p (r = 0) or a 2 (z, r). Shown in Fig. 1 is 
the nonlinear 2D plasma wake wave excited in uniform plasma [N p {r) = 1] by 
the relativistic electron bunch (a — 0.4, a z = 2, a r = 5; for example, in this 
case tibo = 4 x 10 13 cm~ 3 and the characteristic longitudinal and transverse 
sizes of the bunch a z ,r/kp correspondingly are 1.06mm and 2.65mm when 
n p = 10 14 cm -3 ). One can see that the wavelength changes with the radial 
coordinate r. This leads to curving of the phase front and to "oscillations" in 
the transverse direction (see Fig. 2, curve 1). As \z\ increases, the change of 
phase in transverse direction (for fixed z) becomes more and more marked. 
The longitudinal space parameter characterizing the nonlinear wave front 
curving is [10]: 



2[1-A P /A(0)]' 

where A(r) is the nonlinear wavelength. At the distance |Az| « £ from the 
driver the oscillation phase at the axis (r = 0) is opposite to that on the 
periphery (r > o>). Thus, in 2D nonlinear regime the nonlinear wavelength 
changes with r due to nonlinear increase of the wavelength with wave ampli- 
tude. On the other hand, the linear wavelength X p ~ n p decreases with 
equilibrium density of plasma. Hence follows an opportunity to compensate 
the nonlinear increase in wavelength by reducing the wavelength that is due 
to the growth of equilibrium density of plasma. Indeed, assume that the non- 
linear wavelength of the two-dimensional wake wave in the uniform plasma 
A(r) is known. Then, one can roughly compensate for the radial variation of 
the nonlinear wavelength by changing the equilibrium density of plasma in 
the radial direction according to the relation 

A(0)/A(r) = Ap(r)/Ap(0) = K(0)/n p (r)]^. (10) 

In this case the equation for equiphase surfaces is z ~ const, and, therefore, 
the solution for fields could be written in the form fi(z)f 2 (r). If we put that 
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the function A(r) is Gaussian (that is approximately the case at least for 
r < a r , according to numerical data for profiles (8)), then one can take the 
transverse profile of the equilibrium plasma density to be also Gaussian: 

n p (r) = n p0 exp(-r 2 /a 2 p ). (11) 
Then follows from Eqs. (10) and (11) that 

a p = r/MM0)/Mr))] 1/2 . (12) 

For example, according to (12), the numerical data for A(r) in the nonlinear 
wave shown in Fig.l give a p ~ 11. Thus, in the radially nonuniform plasma, 
the density of which is changed according to (11) and (12), one can practically 
avoid undesirable curvature of the wave front of a nonlinear wave. Figs. 3 and 
4 illustrate the validity of this assertion respectively for PWFA and LWFA 
(see also Fig. 2, curve 2). One can see, that the nonlinear wavelength in the 
nonuniform plasma with proper radial profile does not practically change in 
the transverse direction. 

As for the case of LWFA, one has to note that as is well known, without 
optical guiding the diffraction limits the distance of laser-plasma interaction 
(and, hence, the energy gain of particles accelerated by a wake wave) to a few 
Reyleigh lengths Zr = nr^/X, where r the minimum laser pulse spot size at 
the focal point and A is the laser wavelength. For high-intensity laser pulses 
the quantity Zr is usually of the order of several millimeters. One of ap- 
proaches in preventing of diffraction broadening of laser pulse and increasing 
of laser-plasma interaction distance is the guiding of the pulse in preformed 
plasma density channel [3] . Here the unperturbed plasma density grows from 
the pulse axis (r = 0) to its periphery, in contrast with the radial profile of 
the plasma density that was proposed above for prevention of phase front cur- 
vature of a nonlinear wake wave when the plasma density at the driver axis 
is maximum. In case of nonlinear wake wave excitation by high power laser 
pulses (of power P > P c — 2c(e/r e ) 2 [uJ /ujp e (r = 0)] 2 m 17[uj /uj pe (r = 0)] 2 
GW, where r e = e 2 /m e c 2 is the classical electron radius), the process of 
diffraction broadening of a larger part of pulse, including the case of proposed 
equilibrium profile of plasma, may be prevented or essentially retarded due 
to the relativistic self-focusing [3]. One can estimate the condition of rela- 
tivistic self-focusing in our case from the following relation (see, e. g. Sec. 
VI in Ref. [3]): P/P c > 1 + (An/An c )(r s /r ) 4 , where r s is the spot size, 
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An c = \/-Kr e r-Q, An ~ n p (0) — n p (cr r ) is the radial variation of unperturbed 
plasma density in the laser pulse. 

This work has been supported by the International Science and Technol- 
ogy Center under Project No. A-013. 
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FIGURE CAPTIONS 



Fig. 1. The two-dimensional nonlinear wake wave in uniform plasma 
[N p (r) = 1]. The parameters of the bunch are: a = 0.4, a z — 2, a r — 5, 
7 = 10. (a). The density of plasma electrons N e and of the bunch. 1 — the 
density of plasma electrons at the axis, r = 0; 2 — the same for r = 2; 3 - 
r = 4;4 — r = 5;5 — the density of bunch at the axis a(z,r = 0). (b). 
The longitudinal electric field for r = 0, 2, 4 and 5 in the order of magnitude 
reduction, (c). The focusing field f r = (3H e — E r . 1 — r = 2;2 — r = 4;3 
— r = 5. All variables are normalized. 

Fig. 2. The radial behavior of the normalized longitudinal electric field 
strength E z . 1 - - E z (z = — 25, r) in the nonlinear wake wave excited in 
uniform plasma for the case given in Fig. 1 (\Az\ « £ [see Eq. (9)]); 2 - 
E z (z = — 25, r) in nonuniform plasma for the case given in Fig. 3. 

Fig. 3. The two-dimensional nonlinear wake wave in nonuniform plasma 
with a p = 11. The bunch parameters are the same as in Fig. 1. 1 — the 
density of plasma electrons versus z for r = 0,2,4 and 5 in the order of 
magnitude reduction, (b). The same for the longitudinal electric field, (c). 
The focusing field. 1 — r = 2; 2 — r = 4; 3 — r = 5. All variables are 
normalized. 

Fig. 4. The two-dimensional nonlinear wake wave excited by laser pulse. 
The pulse parameters are : Oq = 3.6, a z = 2, a r = 5, 7 = 10. (a). The 
dimensionless accelerating field E z excited in uniform plasma, r = 0,2,4 
and 5 in the order of magnitude reduction, (b). The same in the nonuniform 
plasma, a p = 12. 
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Figure 1 . 
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